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OS | Abstract 

A system of self-gravitating massive fermions is studied in the framework of the 
J> ■ general-relativistic Thomas-Fermi model. We study the properties of the free energy 

functional and its relation to Einstein's field equations. A self-gravitating fermion gas 
£^ | we then describe by a set of Thomas-Fermi type self-consistency equations. 

Thermo dynamical properties of the self-gravitating fermion gas have been extensively 
studied in the framework of the Thomas- Fermi model [1-61. The system was investigated in 

^7s ■ 

— , the nonrelativistic Newtonian limit. The canonical and grand-canonical ensembles for such 
q-( a system have been shown to have a nontrivial thermodynamical limit |l], |2| . Under certain 

conditions this system will undergo a phase transition that is accompanied by a gravitational 
collapse 0, [5| which may have important astrophysical and cosmological implications || [7]. 
>■ . In this paper we formulate the general-relativistic version of the model. The effects 

of general relativity become important if the total rest-mass of the system is close to the 
Oppenheimer-Volkoff limit §. There are three main features that distinguish the relativis- 
tic Thomas-Fermi theory from the Newtonian one: i) the equation of state is relativistic 
ii) the temperature and chemical potential are metric dependent local quantities iii) the 
gravitational potential satisfies Einstein's field equations (instead of Poisson's equation). 

Let us first discuss the general properties of a canonical, self-gravitating relativistic fluid. 
Consider a nonrotating fluid consisting of iV particles in a spherical volume of radius R in 
equilibrium at non-zero temperature. We denote by w M , p, p, n and a the velocity, pres- 
sure, energy density, particle number density and entropy density of the fluid. A canonical 
ensemble is subject to the constraint that the number of particles 

nu^dE^^N (1) 

should be fixed. The spacelike hypersurface £ that contains the fluid is orthogonal to the 
time-translation Killing vector field which is related to the velocity of the fluid 

fc" = e«"; e = (fc%) 1/2 - (2) 
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The metric generated by the mass distribution is static, spherically symmetric and asymp- 
totically flat, i.e. 

ds 2 = £ 2 dt 2 - X 2 dr 2 - r 2 (d6 2 + sin 6d<j) 2 ) . (3) 
£ and A may be represented in terms of the gravitational potential and mass 



t = e*>, (4) 



with 

M(r)= / dr' ^r' 2 p(r') . (6) 
Jo 

The temperature T and chemical potential p are metric dependent local quantities. Their 
space-time dependence may be derived from the equation of hydrostatic equilibrium |J 

d v p=-(p + p)C 1 d v £, (7) 

and the thermodynamic identity (Gibbs-Duhem relation) 

d f = Ud f ~ pd f' ^ 
The condition that the heat flow and diffusion vanish [|K|] 

| = const (9) 

together with (|7|) and (g) implies 

T£ = T ; f4 = fM ) , (10) 

where T and fj,Q are constants equal to the temperature and chemical potential at infinity. 
The temperature T may be chosen arbitrarily as the temperature of the heat-bath. The 
quantity /io in a canonical ensemble is an implicit functional of £ owing to the constraint 
(0). First equation in ( |T0D is the well known Tolman condition for thermal equilibrium in a 
gravitational field [|TT|j . 

Following Gibbons and Hawking [I2| we postulate the free energy of the canonical en- 
semble as 

F = M - J Ta WdEp , (11) 

where M is the total mass as measured from infinity. The entropy density of a relativistic 
fluid may be expressed as 

a = -{p + p- fm). (12) 

Based on equation (|10D the free energy may be written in the form analogous to ordinary 
thermodynamics 

F = M~T S (13) 
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with M = A4(R) and the total entropy S defined as 

S= f R dr47rr 2 \hp + p)-^N, (14) 

Jo 1 1q 

where we have employed the spherical symmetry to replace the proper volume integral as 

/ u"d£„ = ( R drA-nr 2 X. (15) 
js Jo 

The following theorem demonstrates how the extrema of the free energy are related to 
the solutions of Einstein's field equation. 

Theorem 1 Among all momentarily static, spherically symmetric configurations {£(r), A4(r)} 
which for a given temperature To at infinity contain a specified number of particles 

R 

Anr 2 dr X(r)n(r) = N (16) 
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within a spherical volume of a given radius R, those and only those configurations that 
extremize the quantity F defined by ( |13|) will satisfy Einstein's field equation 

d£ = M+ 4nr 3 p 

dr * r(r - 2M) ' 1 ' 

with the boundary condition 

, , / 2M\ 1/2 , x 

Proof. By making use of the identity (§), and the fact that 5(p/T) = 5(pq/Tq) and that N 
is fixed by the constraint fllBD , from equations (|I~3"D and (|14]) we find 

SF = SM- [ R drAnr 2 — (p + p)6X- dr 4nr 2 X—5p . (19) 
Jo T Jo T 

The variations 5X and 5p can be expressed in terms of the variation 5Ai(r) and its derivative 

d5M , 9r 

-- Anr 2 5p. (20) 



dr 

This gives 

5F = 5M- f R dr4nr 2 ^(p + p)^-6M- ["drX^^-. (21) 
Jo 1 oJvl Jo I dr 

By partial integration of the last term and replacing Tq/T by £, we find 



6F=[1- \(R)£(R)] 5M - t dr 

Jo 



SM , (22) 



where 5j\4(r) is an arbitrary variation on the interval [0, R], except for the constraint 
5Ai(0) = 0. Therefore 5F will vanish if and only if 

+ ,)^-£(Afl=0 (23) 
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and 

1 - \(R)£(R) = 0. (24) 

Using (|^) and @, we can write equation (p3|) in the form (fPTD , and equation ( p4|) gives the 
desired boundary condition (|T8[). Thus, 8F = if and only if a configuration {£, M.} satisfies 
equation (|T7D with ([Tj|) as was to be shown. 

Remark 1. A solutions to equation flTT| ) is dynamically stable if the free energy assumes a 
minimum. 

Remark 2. Our Theorem 1 is a finite temperature generalization of the result obtained for 



cold, catalyzed matter [13 



We now proceed to the formulation of the general-relativistic Thomas-Fermi model. Con- 
sider the case of a self-gravitating gas consisting of N fermions with the mass m contained in 
a sphere of radius R. The equation of state may be represented in a parametric form using 
the well known momentum integrals over the Fermi distribution function |T4| 



, d 3 q 1 
71 = 9 Jo (27p 1 + tElT-^IT > ( 25 ) 



roc (flq J7J 

P = 9 L (2nji 1 + e^-,./T • < 26 > 

v = ° T l~W? Hl + e ~ E ' T " IT) > (27) 

where g denotes the spin degeneracy factor, T and fi are local temperature and chemical 
potential, respectively, as defined in equation fllQl) , and E = \Jrri 2 + q 2 . Introducing a single 
parameter 

«=^ = §, (28) 

and the substitution 

£ = ^ , (29) 
m 

equations (p5|) -(^7D may be written in the form 

r°° d 3 g 1 , x 

n H (27T) 3 l + e (^-m)a' ^ 

/-co ^ 
P = H (2^)3 1 + e (E^-m)a ' ( 31 ) 

P = H (27p 3£ 1 + e (^— )« ' (32) 

Field equations are given by 

# _ , + 4?rr 3 p 

rfr _ %(r-2A<) ' [ ' 

dM , 2 



47rr 2 p, (34) 
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with the boundary conditions 



, (R) = ^Mpj' 2 , M ( C)=0. (35) 



Finally, the constraint ([[]) may be written as 



R 







' dr 47rr 2 (l - 2M /r)" 1/2 n{r) = N. (36) 



Given the ratio a, the radius R, and the number of fermions N, the set of self-consistency 
equations (|30|)-(|3"6"|) defines the Thomas- Fermi equation. One additional important require- 
ment is that a solution of the self-consistency equations (|30l) -(|36"l) should minimize the free 
energy defined by (|D|). 

We now show that a solution of the Thomas-Fermi equation exists provided the number 
of fermions is smaller than a certain number iV max that depends on a and R. From ( |31~1) 
and (|32|) it follows that for any a > 0, the equation of state p(p) is an infinitely smooth 
function and dp /dp > for p > 0. Then, as shown by Rendall and Schmidt fl5|], there exist 



for any value of the central density po a unique static, spherically symmetric solution of field 
equations with p — > as r tends to infinity. In that limit Ai(r) — > oo, as may easily be seen 
by analysing the r — > oo limit of equations (|3"3] ) and fl3"4"|). However, the enclosed mass M 
and the number of fermions N within a given radius R will be finite. We can then cut off 
the matter from R to infinity and join on the empty space Schwarzschild solution by making 



use of equation (|35|). This equation together with fl28|) fixes the chemical potential and the 
temperature at infinity. Furthermore, it may be shown that our equation of state obeys a 
7-low asymptotically at high densities, i.e., p = const n 7 and p = (7 — l)p, with 7 = 4/3. It 
is well known |13| that in this case, there exist a limiting configuration {ipoo{ r ), -M( r )oo} such 
that M and iV approach non-zero values and iVoo, respectively, as the central density 
po tends to infinity. Thus, the quantity iV is a continuous function of po on the interval 
< Po < °°> with iV = for p = 0, and iV = as p — > 00. The range of iV depends 
on a and R and its upper bound may be denoted by N m ^(R, a). Thus, for given a, R and 
iV < N max (R, a) the set of self-consistency equations nas a t least one solution. 

Next we show that, in the Newtonian limit, we recover the nonrelativistic Thomas- Fermi 
equation. Using the nonrelativistic chemical potential p^R = Po — m and the approximation 
£ = e </3 ~l + <p, -E~m + q 2 /2m and Ai/r 1 , we find the usual Thomas-Fermi self- 
consistency equations |5|, |j 

p r° d 3 q ( q 2 m Pnr.Y 1 

n = — = g —— 3 1 + exp — — + — <p - ^— , 37 
m Jo [2TT) d \ 2mT T T J 

dM , 2 



d 3 q 


(■ 


(2tt) 3 




dip 




dr 




<p(R)-- 





dr 

mN 



R 

R 

2 



4vrr>, (38) 
M(0) = 0, (39) 



drAirr n(r) = N. (40) 
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The free energy (|T3| ) in the Newtonian limit yields 



1 r R f R 

F = mN + unrN — — / dr Anr 2 mp - / dr4irrp (41) 

2 Jo Jo 

with 

[°° d 3 q f q 2 m Hnrs\ (ao , 

which, up to a constant, equals the Thomas- Fermi free energy ||. 

A straightforward thermodynamic limit iV — > oo as discussed by Hertel, Thirring and 
Narnhofer || |J is in our case not directly applicable. First, in contrast to the non-relativistic 
case, there exists, as we have demonstrated, a limiting configuration with maximal M and N. 
Second, the scaling properties of the relativistic Thomas-Fermi equation are quite distinct 
from the nonrelativistic one. The following scaling property can be easily shown: If the 
configuration {ip(r), J\4(r)} is a solution of the self consistency equations (|30|)-(j36|), then the 
configuration {ip = ^(A^r), M. = AM. (A^r); A > 0} is also a solution with the rescaled 
fermion number iV = A 3 ^ 2 N, radius R = AR, asymptotic temperature T = A _1//2 T , and 
fermion mass rh = A~ l / 2 m. The free energy is then rescaled as F = AF. Therefore, there 
exist a thermodynamic limit of N~ 2 ^ 3 F, with N~ 2 ^ 3 R, A^ 1//3 T , N 1 ^ 3 m approaching constant 
values when N oo. In that limit the Thomas- Fermi equation becomes exact. 

It is obvious that application of this model to astrophysical systems should work very 
well if the interactions among individual particles are negligible. This applies, for example, 
to weakly interacting quasidegenerate heavy neutrino or neutralino matter [6,7,16-19]. or 



perhaps even to collisionless stellar systems |2Q, 21 
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